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Abstract 

A new mechanism of cosmic axion production is proposed: axion bremsstrahlung from 
collisions of straight global strings. This effect is of the second order in the axion coupling 
constant, but the resulting cosmological estimate is likely to be of the same order as that 
corresponding to radiation from oscillating string loops. This may lead to a further restriction 
on the axion window. 



1 Introduction 

The cosmic axion remains one of the viable candidates for dark matter, though the range 
of masses, which remains open accounting for both collider and cosmological constraints, is 
rather narrow. The model contains an unknown parame ter , th e vacuum expectation value /, 
marking the energy scale of the U(l) symmetry breaking ' 1 ^ 1 The mass of the axion, which 
is acquired after the QCD phase transition, is inversely proportional to /, whose upper bound 
is of cosmological origin and follows from the requirement that axions produced during the 
cosmological evolution do not overdose the Universe. The lower bound on the axion mass 
lies between several units and several tens of fieV ; this corresponds to the value of / between 
10 12 and 10 11 GeV. For a rec ent survey of the present theoretical and astrophysical status 
of the cosmic axion model se e^^t earlier reviews included 

The axion string networkl^l^l is formed at the temperature of the Peccei-Quinn phase 
transition /, and it is usually assumed that the reheating temperature is higher than /, 
otherwise the network would be diluted by inflation. Strings are primarily produced as 
long straight segments whose length i s of the order of the horizon size, and they initially 
move with substantial friction^ElE] d ue t the scattering of cosmic plasma particles. At 
some temperature T* < / scattering becomes negligible, and the string network enters the 



scale-invariant regi me P^ ^^ when strings form closed loops and move almost freely with 
relativistic velocitie d^^ . The standard estimates of the axion radiation from global strings 
are based on the assumption that the main contribution comes from the oscillating string 
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Here we discuss another mechanism for axion radiation: bremsstrahlung, which has to be 
produced in collisions of long strings. The existence of such an effect can be demonstrated as 
follows. Consider two infinite straight strings inclined with respect to each other and moving 
in parallel planes. Due to interaction via the axion field, strings will be deformed around 
the point of minimal separation between them. The motion of this point is not associated 
with the propagation of any signal, and the corresponding velocity may be arbitrary, in 
particular, superluminal. In the second order of string-axion interaction the superluminally 
moving deformation must produce Cerenkov axion radiation. A similar mechanism was 
earlier suggested for gravitational radiation of local strings, but in that case the explicit 
calculations'^ have led to a zero result. The vanishing of gravitational radiation, however, 
has a specific origin related to the absence of gravitons in 2 + 1 gravity theory. Indeed, 
as was explained inEB, two crossed superluminal strings can be "parallelized" by suitable 
coordinate transformations and world sheet reparamctcrizations, so that the collision of 
strings is essentially equivalent to the collision of point particles in 2 + 1 dimensions. In the 
case of the global strings similar considerations lead us to the problem of the electromagnetic 
radiation of point charges in 2 + 1 dimensions, which is not forbidden by dimensionality. Thus 
one can expect a non- vanishing axion bremsstrahlung from collisions of global strings. 

2 Strings interacting via axion field. Axion radiation. 

Consider a pair of relativistic strings x^ = x^ l (a a ) 1 /Lt = 0,1,2,3, a a = (r, a), a = 0, 1, 
where n = 1,2 is the index labelling them. The 4-dimensional space-time is assumed to be 

flat and the signature is +, (and (+, — ) for the string world-sheets). Strings interact 

via the axion field B^ v as described by the action^ 



S = - >^ / {^V^l ab daX^d b <V^ + ^fnB^e ab d a x^d b x v n \d 2 a. 




+ \ / H^ x H^ x d 4 x. ! I ) 

Here are the (bare) string tension parameters, /„ are their axion couplings (n-labelling 
helps to control the perturbation expansions), the Levy-Civita symbol is chosen as e 01 = 
1, "f a b is the metric on the world-sheets. The totally antisymmetric axion field strength 
is defined as = d^B u \ + d u B\^ + d\B^ u with the Lorentz gauge d\B vX = 0. The 

corresponding potential two-form is the sum B^ v = B^ v + Bt^ of contributions B£ v due to 
each string: 

rf p d a d p B^ = -47rJr, (2) 

where J%» = If J V£ v 5\x - x n {a n )) d 2 a n and V£ v = e ab d a x£d b x" n . 

In the gauge j a b = Vab the renormalized strings equations of motion read : 

M „(^-^)<=27r/ n ff^V i r A ' ) n^n', (3) 

where the field corresponds to contribution of the n'—th string (no sum over n'). The 

corresponding constraint equations are x 2 + x' 2 = 0, x^x^rj^ = for each string, where 
dots and primes denote the derivatives over r and a. The self-action terms in the equations of 
motion diverge both near the strings and at large distances, so two regularization parameters 
5 and A have to be introduced!^ which are absorbed by the classical renormalization of the 
string tension as 

/i = /i + 2^/ 2 log(A/(5). (4) 



The ultraviolet cutoff length 5 is of the order of the string thickness S ~ 1//, while the 
infrared cutoff distance A usually is chosen as the correlation length of the string network. 



Assuming that such a renormalization is performed, we are left with the same equations 
of motion |JSJ, now with the physical tension parameters in which only the mutual 
interaction terms should be used in the right hand side. 

Our calculation follows the approach of^J and it consists in constructing solutions of 
the string equations of motion and the axion field iteratively in terms of the string-axion 
field interaction constant (equal to /). We would like to mention also a similar perturba tive 
approach to solve the problem of gravitating point masses in terms of geodesic deviations^!. 
In the zero order approximation the strings are moving freely, so the lowest order axion fields 
due to both strings describe their mutual interaction and do not contain the radiative 

part. Substituting them to we obtain the first order deformations of the world-sheets. 
These are used to build the first order source terms in the axion field equations J2J) 

which generate the second order axion fields already containing the radiative parts. 

The radiation power can be computed as the reaction work produced by the half sum of the 
retarded and advanced fields upon the source and presented in the standard form^ 

E = \ J J aP {k)\*5{k 2 )d% (5) 

where e is the sign function. 

The final formula for the axion bremsstrahlung from the collision of two global strings can 
be obtained analytically in the case of the ultrarelativistic collision with the Lorentz factor 
7 = (1 — t> 2 ) -1 / 2 3> 1. The spectrum has an infrared divergence due to the logarithmic 
dependence of the string interaction potential on distance, so a cutoff length A has to be 
introduced. The final formula for the radiated energy per unit length of the target string 
reads: 

d£ 167r 5 /V d 

~dl = ~ 3~2 — y ~ ^Ja' K = ^ cosa ' W 
here a is the strings inclination angle and 

*<»>=VMiHH- 31n(%eC)+ ? 

where the generalized hypergeometric function is introduced and C is the Euler's constant. 
The leading term for small y is the logarithm, while the asymptotic value for large y is 
constant </>(oo) = 7/2. 

The velocity of the effective source (i.e. of the point of minimal separation between the 
strings) tends to infinity when the string become parallel. In this case, corresponding to 
k = 7 in the above formula, the Cerenkov cone opens up to 9 — tt/2, and the whole picture 
becomes essentially 2 + 1 dimensional. This gives an alternative description of the effect as 
bremsstrahlung under collision of point electric charges in 2 + 1 electrodynamics^ 1 . 



3 Cosmological estimates 

Our mechanism presumably works in the temperature interval (T*, Ti) between the onset 
of the scaling regime till the QCD phase transition, though an additional contribution form 
the damped epoch T < T < can also be expected. The scaling density of strings is 
determined from numerical experiments, it can be presented as p s — (,[i/t 2 with £ varying 
from 1 to 13 in different simulations. The corresponding moments of time are commonly 
estimated as^l 

10-20 (loifev) ^ <l(sec ^ 2 ' 10-7 (i^4v) 1/3 - (7) 

Consider the scattering of an ensemble of randomly oriented straight strings on a selected 
target string in the rest frame of the latter. Since the dependence of the string bremsstrahlung 
on the inclination angle a is smooth, we can use for a rough estimate the result for parallel 
strings (a — 0) introducing an effective fraction v of "almost" parallel strings (roughly 1/3). 
Then, if there are N strings in the normalization cube V = L 3 , we have to integrate the 



radiation energy released in the collision with the impact parameter d — x over the plane 
perpendicular to the target string using the measure N/L 2 ■ 2irxdx. Actually we need the 
radiation power per unit time, so, for an estimate, we have to divide the integrand on the 
impact parameter. This quantity should be multiplied by the total number of strings N to 
get the radiation energy released per unit time within the normalization volume. Therefore, 
for the axion energy density generated per unit time we obtain: 



de a _ f L N_N 

dt ~ L tv v v 



Sv-^— 2-ixdx. (8) 



Substituting here the Eq. JSJ, we integrate over the impact parameter, taking into account 
that the string number density is related to the energy density via N/V — p s j fj,L. Then 
introducing the integral £(w) = J <p(y)dy and keeping only the second (leading) term in the 
expression for the string tension fi, we obtain 

where K = |-7r 4 r/£ 2 . Since in the cosmological context L ~ A ~ t, one can take w = j~ 2 
for an estimate. The "realistic" value of 7 is of the order of unity, while our formulas were 
obtained in the 7> 1 approximation. But an independent calculation shows that for 7 ~ 1 
the contribution from the first term is of the same order, so we can hope to get reasonable 
order of magnitude estimate using the above formulas for 7 not very large. The full 7- 
dependence is given by the function 7 3 £(7 -2 ), which for 7 = \[2 is 7 3 £(7~ 2 ) ~ 8.137, while 
for 7 = 5,7 3 £(7~ 2 ) « 55.83. Another source of uncertainty is the numerical value of the 
parameter £, which was obtained in different simulations within the region between 1 and 
13. 

Integrating the radiation power over time from t* to t\ as given by Q one finds the 
energy density of axions at the moment of QCD phase transition t\ ^> t*: 



2tl ]n 2 (U/to)' 



(10) 
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Finally dividing this by the critical energy density at t = t\ e cr = 32 Jt[ , wipi = 1.22 • 
10 19 GeV, and using the relation <*/t = (rnpi//) 2 , we obtain for 7 = y/2 the following 
estimate for the relative contribution of the bremsstrahlung axions 

ft!f - 0.5 10 16 (£\ ( „L ' • (11) 



13/ V 10l2GeV 

This rough estimate shows that our new mechanism is very sensitive to the value of / and 
gives an upper bound on the axion window of the order of several 10 10 GeV. We postpone 
more accurate cosmological estimates for a separate publication. In particular, more careful 
analysis is needed to accommodate the analysis to the actual astrophysical value of the string 
Lorentz factor. 
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